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Order parameter fluctuations (the largest cluster size distribution) are studied within a three- 
dimensional bond percolation model on small lattices. Cumulant ratios measuring the fluctuations 
exhibit distinct features near the percolation transition (pseudocritical point), providing a method 
for its identification. The location of the critical point in the continuous limit can be estimated 
without variation of the system size. This method is remarkably insensitive to finite-size effects and 
may be applied even for a very small system. The possibility of using various measurable quantities 
for sorting events makes the procedure useful in studying clusterization phenomena, in particular 
nuclear multifragmentation. Finite-size scaling and A-scaling relations are examined. The model 
shows inconsistency with some of the A-scaling expectations. The role of surface effects is evaluated 
by comparing results for free and periodic boundary conditions. 
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I. INTRODUCTION 

The main motivation of nuclear multifragmentation 
studies is probing a liquid-gas coexistence region in the 
phase diagram of nuclear matter Q,H0| . Many works de- 
duce the occurrence of a first- or second-order phase tran- 
sition iisnsssHmmmH. Aith ° u g h b ° th 

transition types can be expected, unambiguous identifica- 
tion is difficult due to strong finite-size effects in systems 
with a small number of constituents. In such systems, for 
example, a first-order phase transition may mimic critical 
behavior [15lll6l| . On the other hand, nuclear multifrag- 
mentation induced by high energy collisions shows strik- 
ing similarities to percolation processes which are known 
to contain a second-order phase transition (critical be- 
havior) mHHHIIlEllillH. Percolation - based 

models seem to be also successful in describing fragmen- 
tation of atomic clusters [24l l25j . These observations 
lead to formulation of the hypothesis that percolation 
could be a universal fragmentation mechanism for simple 
fluids [23]. For better recognizing critical-like behavior 
observed in fragmentation processes, simultaneous appli- 
cation of various complementary methods is necessary. 
Percolation models are often used to construct or verify 
procedures tracin g cr itical behavior in fragmenting sys- 
tems [H H E 113, IM HI HI |H. They provide a 
simple tool for studying universal aspects of the critical 
behavior and the role of finite-size effects. 

In the framework of a percolation model we have ex- 
amined the largest cluster size distribution. The size of 
the largest cluster, as an order parameter in aggrega- 
tion scenarios of the fragment production, is of partic- 
ular interest in phase transition studies [28l l29j| . The 
limiting forms of the distribution for normal phases are 
predicted by classical limit theorems for random vari- 
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ables [H |33, I^H \3^. At a second-order phase tran- 
sition the system is highly correlated with fluctuations 
occurring on all length scales. Properties of the or- 
der parameter close to the critical point can be stud- 
ied with the renormalization group and finite-size scal- 
ing approaches [23, [23, LLJ L^s uM- Botet and 
Ploszajczak have proposed to identify the second-order 
critical behavior in finite systems by examining univer- 
sal features of the order parameter fluctuations with a 
A-scaling method H |H El H3 ■ The method has been 
ap plie d to several models and nuclear fragmentation data 
B Hi GB III [11 III III III II| • We will confront A- 
scaling predictions with percolation results. 

In order to compare theoretical predictions with frag- 
mentation data all experimental conditions should be 
carefully considered. The bulk behavior of the order pa- 
rameter can be significantly modified in small fragment- 
ing systems by finite-size and boundary effects. The con- 
trol parameter is usually not well measured and must 
be substituted by other measurable quantities in sort- 
ing events, leading to additional modifications. We aim 
to evaluate the significance of such effects in the present 
work. 

The calculations were performed with a three- 
dimensional bond percolation model on the simple cubic 
lattices |33j . A Monte Carlo procedure with the Hoshen- 
Kopelman cluster labeling algorithm was employed to 
generate events for a distribution of the bond probabil- 
ity, p, being the control parameter. The lattices of size 
N = L 3 with L = 3, 4, 5, 6 correspond to the range of sys- 
tem sizes available in nuclear reactions. Free boundary 
conditions were applied to account for the presence of sur- 
face in real systems. To evaluate the role and importance 
of finite-size and boundary effects the calculations were 
extended to include larger systems and periodic bound- 
ary conditions. 

We study low-order cumulants (cumulant ratios) as the 
mean, variance, skewness and kurtosis of the largest clus- 
ter size distribution. These standard statistical measures 
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contain the most significant information, providing a ro- 
bust identification of the percolation transition. In order 
to place our results in a wider context, we will briefly 
recall in the next section some signals of criticality that 
are frequently tested in fragmentation studies. 



II. PERCOLATION TRANSITION IN SMALL 
SYSTEMS 

In the strictest sense, the phase transition occurs in 
the continuous limit N — > oo. Then, below the percola- 
tion threshold p < p c ~ 0.2488, only finite clusters are 
present. When p > p c there exists an infinite cluster 
spanning the whole lattice. The fraction of sites belong- 
ing to this cluster is the order parameter. In finite sys- 
tems the transition is smoothed. The probability that at 
least one cluster connects the bottom and the top lattice 
planes changes gradually as illustrated in Fig. 1(a). A 
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FIG. 1: Bond percolation on a lattice of linear size L = 5 with 
free boundary conditions. Plotted as a function of the bond 
probability: (a) the probability that a lattice is spanned, (b) 
moments of the average distribution of cluster sizes, (c) the 
variance of the largest cluster size distribution with various 
normalizations, (d) The average distribution of cluster sizes. 

natural way to locate the transition is to examine quan- 
tities which diverge at the critical point in the continu- 
ous limit. For finite systems, the divergence is replaced 
by a maximum located near p c . This can be seen in 
Fig. 1(b) for the second moment of the average distribu- 
tion of cluster sizes, M2, the mean cluster s ize, M2/M1, 
and the reduced variance, 72 = M 2 M /M 1 2 j20|. M k de- 
notes the fc-th moment of the cluster size distribution, 
Mk = Y] s s k n s , where n s is the average number of clus- 



ters of size s normalized to the system size (the largest 
cluster excluded). The mean cluster size is the analog of 
the susceptibility and the location of its maximum value 
defines the pseudocritical point. Another frequently used 
method is a power-law fit to the fragment size distribu- 
tion 0, [22I ■ In our example the best fit is obtained for 
p = 0.33 as shown in Fig. 1(d). The fragment size dis- 
tribution follows in some range the asymptotic behav- 
ior n s = qos~ T , where the Fisher exponent r = 2.189 
|39| . The normalization constant qo = 1/ •s^ 1 ~ r - ) 
is taken as 0.173 from the summation computed up to 
s = 10 9 . Another example is the maximum fluctuation 
of the largest cluster size. The size of the largest clus- 
ter, s max , plays the role of an extensive order parameter. 
Its fluctuations are usually measured by the variance, K2, 
or the normalized variance, ty-il ( sm. q x), of the probability 
distribution P(s max ) [2(1 [23, 123, E3 The two quantities 
are peaked in the transition region as shown in Fig. 1(c) 
by the dotted and dashed lines. When using the normal- 
ization K2/ {smax) 2 , the maximum is located remarkably 
close to the critical point. 

The above examples show that various investigated sig- 
nals appear at different positions and in most cases they 
are shifted from the critical point toward the ordered 
phase region. In small systems the shifts may be signif- 
icant and should be taken into account in a criticality 
analysis. The location of the critical region in nuclear 
multifragmentation is often deduced from a power-law 
fit to the fragment size distribution. This location, ap- 
pearing near the pseudocritical point, would correspond 
to a temperature T pc distinctly lower than the true crit- 
ical temperature T c . For example, converting the bond 
probability to the temperature with the prescription of 
Ref. [HEl one obtains T pc /T c = 0.64,0.73 and 0.78 
for N = 27,64 and 125, respectively. 

The location of the true critical point, p c , is of par- 
ticular interest. According to the finite-size scaling the 
position of a signal, p(L), is expected to converge to p c 
with increasing linear lattice size L as 

pi^-p^L- 1 '", (1) 



where v is the correlation length exponent. Estimation 
of a critical point by such an extrapolation method seems 
to be difficult in the case of nuclear multifragmentation. 
Sizes of fragmenting systems created in nuclear reactions 
are not well controlled due to the preequilibrium emission 
and their range is limited. In addition, one may expect 
large departures from the scaling relation of Eq. (1) for 
very small systems. Our observation is that, without re- 
lying on the finite-size scaling, the best estimation of the 
critical point is given by the position of the maximum of 
K2/ (smax) 2 ■ Behavior of this quantity for different sys- 
tem sizes and various event sortings will be investigated 
in the following sections. 
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III. ORDER PARAMETER FLUCTUATIONS 
A. Cumulants and finite-size scaling 

The order parameter probability distribution represen- 
tative for small lattices with open boundaries is shown in 
Fig. 2 for various values of p. Far from the transition 
the distribution is sharply peaked with an extended tail 
to the right (left) in the disordered (ordered) phase and 
positioned close to the limiting values. In the transition 
region the distribution rapidly evolves passing through a 
broad, flattened and (nearly) symmetrical distribution. 
This behavior can be well characterized by using the 
skewness, K 3 , measuring the asymmetry of a distribution 
and the kurtosis excess, K±, which quantifies the degree 
of peakedness. The quantities of interest are defined and 
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FIG. 2: Probability distributions of the largest cluster size for 
a lattice of size N = 125 with free boundary conditions. 



expressed in terms of the cumulants as 

K 2 = /i 2 /(s mal ) 2 =k 2 /k\ 



M3/M2 /2 

m/i4 - c 



I 3/2 
= K3/K2 



K4/K. 



(2) 



where m — {{s max — {s max )) 1 ) is the i-th central moment, 
and Hi is the i-th cumulant of the P{s max ) distribution. 
They are plotted as a function of p for different system 
sizes in Fig. 3(a). In the vicinity of the critical point and 
for L — > 00 one expects for these dimensionless parame- 
ters the scaling relation 



Ki = fi[(p-p c )L 



(3) 



If the scaling holds, values of Ki at p c are independent of 
the system size, which is approximately observed in our 
plots as the crossing of curves for different L near p c . To 
verify the scaling in the neighborhood of the critical point 
the Ki parameters are replotted in Fig. 3(b) against 
the scaling variable (p —p c )L x / v , where v = 7/8. The 
collapse of the data shows that the scaling relation with 
no corrections for finite-size effects is well satisfied even 




FIG. 3: The cumulant ratios of Eq. (2) as a function of the 
bond probability (a) and the scaling variable (b) for systems 
of different sizes. Calculations with free boundary conditions. 



for such small lattices with open boundaries. As can be 
estimated from Fig. 3(b) the asymptotic values of Ki at 
p c are about K 2 = 0.24, K 3 = 1.2, and K4 = 1.5. For 
small systems they are somewhat smaller with largest 
deviations observed for K 2 . A prominent feature of the 
K 2 distribution is the maximum located very close to p c 
irrespective of the system size. 

Another characteristic point, corresponding to the 
broad transitional distribution shown in Fig. 2, is where 
^3 = and K4 reaches its minimum value of about —0.9, 
and K 2 — 0.135. It is observed at some distance from 
the critical point, which depends on the system size as 
(p ~ Pc) — O-SL^ 1 /". This point approximately coin- 
cides with the maximum of the mean cluster size and the 
power-law behavior of the fragment size distribution (see 
Fig. 4), and may be used as an estimation of the pseudo- 
critical point. The line in Fig. 4 is a power-law fit of Eq. 
(1) to such points giving v = 0.878 ± 0.005 in agreement 
with the expected value v = 0.875. Corrections to the 
scaling are not significant in this case. Other variables 
considered in Fig. 4 show much larger deviations from 
the asymptotic scaling behavior. 

The above characteristics are for free boundary condi- 
tions. Behavior of Ki for periodic boundary conditions 
when surface effects are reduced is shown in Fig. 5. Also 
in this case the finite-size scaling features are clearly ob- 
served. As expected, the scaling functions fa are com- 
pressed now toward lower bond probabilities, which can 
be seen by comparing Fig. 5(b) with Fig. 3(b). The 
pseudocritical point is positioned very close to the criti- 
cal point. Thus, the large difference in locations of these 
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points in the case of free boundary conditions may be 
interpreted as a surface effect. 
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FIG. 4: Finite-size scaling plot for values of the bond proba- 
bility, p, at which the conditions indicated on the figure are 
fulfilled. The line is the best linear fit to the open circles. 
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FIG. 5: Same as Fig. 3 for periodic boundary conditions. 
The dots in the left part indicate the Kz and K4 values for 
the Gumbel distribution (see text). 



B. Delta scaling 

Behavior of the cumulant moments is of interest in the 
context of A-scaling prop osed for studying criticality in 
finite systems 0, Ha, I29L l3fjj . Probability distributions 




FIG. 6: The normalized variance of P(s max ) versus the bond 
probability for free boundary conditions (a), and periodic 
boundary conditions (b). 



P(smax) of the extensive order parameter, s maxi for dif- 
ferent "system sizes", (s max ), obey A-scaling if they can 
be converted to a single scaling function $(z(m) by the 
transformation 



{Smax) A P{s m ax) = $0(A)) = * 



\^max / 



(4) 

where 1/2<A<1. ft was argued that the order param- 
eter satisfies the "first-scaling law" (A = 1 scaling) for 
systems at the critical point and in the disordered phase. 
The "second-scaling law" (A = 1/2) applies for systems 
in the ordered phase far from the transition point. Us- 
ing a three-dimensional bond percolation model, it was 
demonstrated that the A = 1 scaling holds near the crit- 
ical point while the A — 1/2 scaling is satisfied above the 
percolation threshold atp = 0.35 |23,|2^|. However, these 
results were obtained for rather large systems N = 14 3 
to 32 3 with periodic boundary conditions. In the follow- 
ing, we examine the scaling properties in a wide range 
of the control parameter, also for smaller systems with 
open boundaries. 

In case of a A-scaling the normalized cumulants Ki / k\ a 
are independent of the "system size", K\ = (s max ). 
Therefore, for a set of P(s max ) distributions, K3 — const 
and = const are necessary conditions for a A-scaling. 
Additionally, K2 — const for A = 1, and K2/K1 = const 
for A = 1/2. The conditions for the A = 1 scaling are 
fulfilled only in the vicinity of the critical point. The 
crossing points seen in Figs. 3(a) and 5(a) indicate the 
system size independence of Ki near p c . The conditions 
are also satisfied around the critical point when both the 
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system size and the control parameter are varied so that 
(p — PcjL 1 '" = const, accordingly to the finite-size scal- 
ing. Since the conditions are necessary but not sufficient, 
we have checked that indeed the scaling relation of Eq. 
(4) is approximately satisfied. Observing the A = 1 scal- 
ing requires then a variation of the system size without 
or with a very specific change of the control parameter. 
It cannot be observed when the system size is fixed. A 
similar conclusion has been reached in Ref. (3^ for a lat- 
tice gas model. This contradicts the statement of Refs. 
[28l 123, l30| that the scaling relation is valid indepen- 
dently of any phcnomcnological reasons for changing the 
"system size" . 

Fig. 5(a) shows no evidence for the presence of a 
A-scaling in the subcritical region p < p c (disordered 
phase) . The largest cluster size in subcritical percolation 
have been extensively studied in Ref. j^H. As predicted 
by the theory of extremes of independent random vari- 
ables, P(s m ax) converges to the Fisher-Tippett (Gum- 
bel) distribution when N — > 00. The mean grows log- 
arithmically with the system size while the variance is 
bounded. Such a behavior cannot be described by a A- 
scaling. For the Gumbel distribution K3 ~ 1.14 and 
K4 = 2.4, marked in Fig. 5(a) by the dots. As can be 
seen, the small systems show significant deviations from 
these asymptotic values even for periodic boundary con- 
ditions. 

The limiting behavior of the largest cluster size in the 
supercritical region p > p c (ordered phase) is governed 
by the Central Limit Theorem [29l l32| . The asymptotic 
distribution is Gaussian, K3 = K4 = 0, with the mean 
and variance both increasing linearly with the system 
size, K2/K1 — const, satisfying the A = 1/2 scaling rela- 
tion. Such characteristics are seen in Figs. 5(a) and 6(b) 
away from the critical point for larger systems with peri- 
odic boundary conditions. Since the normalized variance 
K2/K1, same for all N, systematically decreases with p, 
the scaling can only be observed for fixed p and differ- 
ent N. Considering very small systems L — 3 to 6 with 
free boundary conditions as appropriate for nuclear ap- 
plications, Fig. 6(a) shows that K2/K1 a t fixed p strongly 
depends on the system size. This indicates the violation 
of the scaling as a consequence of surface effects. 

Investigations of the largest fragment charge distri- 
bution, P(Z max ), observed in heavy-ion central colli- 
sions at bombarding energies between 25 and 150 AMeV 
have shown that K^/ni — const at lower energies while 
K2 — const in a high energy range 0, H^, |3(J . The two 
regimes appear on the ln(/t2) versus \n.{n\) plot along two 
lines with the slope of 1/2 and 1. This observation has 
been interpreted as the presence of the two limiting A- 
scaling laws corresponding to the ordered and disordered 
phases (in the case of A-scaling the slope is equal to 
A). Fig. 7 shows such a plot when experimental events 
are sorted according to the estimated source excitation 
energy. These data taken from Ref. [3!| include quasi- 
projectiles from Ta + Au collisions allowing to observe 
a strong suppression of the fluctuations at lowest ener- 
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FIG. 7: Log-log plot of the variance versus the squared mean 
value of the largest fragment size (charge) distribution. The 
points are experimental data for events sorted according to 
the excitation energy [3f|. The lines are percolation results 
for TV = 64 when events are sorted by the bond probability 
(dashed line) and by the fraction of open bonds (dotted line). 



gies. Assuming the percolation pattern of the cumulants, 
interpretation of the high energy branch in terms of the 
A = 1 scaling would require that for different excita- 
tion energies all fragmenting systems are created with 
nearly the same value of a control parameter close to a 
critical point. More realistically, the control parameter 
varies with the excitation energy while changes of the 
system size are less significant. The dashed and dotted 
lines in Fig. 7 show percolation results for fixed N = 64 
when events are binned by the bond probability and by 
the fraction of open bonds. Given the mean, the vari- 
ance depends not only on the system size but also on the 
choice of binning variable. For a quantitative comparison 
with the data one would have to determine appropriate 
system sizes and a sorting variable equivalent to the ex- 
citation energy. Nevertheless, the qualitative behavior 
of the lines shows similarity to the experimental data. 
The model suggests that the slope changes continuously 
and, as is shown in Fig. 8, the point with slope of 1 cor- 
responds to the maximum of K2 (locally K2 — const), 
whereas the slope of 1/2 reflects the maximum of K2/ n\- 
These features are not related to a A-scaling. The rise 
and fall behavior of the correlation in Fig. 7 is a simple 
consequence of the mass conservation constraint. Some 
points on this line may have a particular meaning de- 
pending on the assumed model. In the present model, 
the point of slope 1 approximately corresponds to the 
critical point. Within the canonical lattice gas model the 
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FIG. 8: Percolation lattice of 64 sites with open boundaries, 
(a) The normalized variance of P(s max ) as a function of the 
bond probability, (b) The variance versus the squared mean 
in log-log representation for events binned by p. 



point of maximum variance occurs close to the critical 
point at the critical density. For subcritical densities it 
is located inside the coexistence region [38) . 
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FIG. 9: The cumulant ratios as a function of the fraction of 
open bonds; on top the correspondence to the mean value of 
the bond probability. 



C. Event sorting effects 

The cumulant properties of the largest cluster size 
distribution presented above are for an ideal situation, 
which assumes that generated events are sorted accord- 
ing to precisely known values of the control parameter. In 
experimental studies such a selection is difficult to realize. 
Usually the sorting parameter is a measurable quantity 
(e.g., multiplicity, excitation energy per nucleon) which 
is correlated with the control parameter with a signif- 
icant dispersion. Even if an attempt is made to esti- 
mate a control parameter such as the temperature in 
nuclear multifragmentation, some dispersion is unavoid- 
able. This can be simulated by using the fraction of open 
bonds, q, for sorting events. On average, the relation be- 
tween p and q is linear, q ~ 1 — p, for all system sizes. 
Fig. 9 presents Ki as a function of q for small lattices 
of L = 3 to 6 with open boundaries which are relevant 
for nuclear multifragmentation studies. The top diagram 
shows the correspondence between q and (p) . Comparing 
Ki(q) with Ki(p) in corresponding intervals shows that 
absolute values may change significantly, however, some 
characteristic features are approximately preserved. The 
2^2 (?) distributions exhibit maxima near the "critical" 
value q c ~ 0.75 corresponding to p c . The zeros of K%(q) 



coincide with the minima of K4 (q) reflecting the behavior 
near the pseudocritical point. However, at the "critical" 
point, q Cl the cumulant values are now much smaller and 
differ with the lattice size. The crossing points appear 
shifted from q c toward the ordered phase. These inter- 
section points are spread out over some range of q and 
are expected to converge to q c with increasing N. The 
relation between p and the number of broken bonds is 
governed by a binomial distribution. This implies that 
for a given q, the dispersion of p vanishes as ~ ./V" 1 / 2 with 
N — > 00. On the other hand, according to the finite-size 
scaling, Ki(p) follows the same pattern irrespective of N 
within a fixed interval of the scaling variable (p—p c )L 1 / v . 
Thus, the corresponding interval of p vanishes with in- 
creasing TV as ~ N~ 1 / 31 ' ~ tv -1 / 2 - 625 . Since the disper- 
sion vanishes faster, the limiting distributions Ki(q) and 
Ki(p) will be equivalent: Ki(q) — > Ki(p = 1 — q) when 
N — > co. Significant differences between the sortings are 
observed in small systems. 

Particularly interesting is sorting events according to 
directly measurable quantities. We have examined the 
Ki dependencies on the overall multiplicity, m, (Fig. 10) 
and the total size of all clusters of size greater than 1, 
Sbound, (Fig. 11). The parameters are normalized to 
the system size N . In both cases the average correspon- 
dence to the control parameter depends on the system 
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FIG. 10: The cumulant ratios as a function of the normalized 
multiplicity. 



FIG. 11: The cumulant ratios as a function of the normalized 
total size of complex fragments. 



size, which is shown in the top diagrams where the solid 
line is for L = 6 and the dotted one for L = 3. The 
K2 distributions exhibit maxima whose positions corre- 
spond to the critical value of the control parameter, p c . 
For such "near-critical" events, values of the sorting vari- 
ables and the cumulant ratios depend on the system size 
as shown in Fig. 12. They all can be well described by 
the equation 



with coefficients a, 6, and c listed in Table I. In practice, 
when the system size in not well known, one can examine 
Ki as a function of m (and/or Sbound)- The system size 
N can be found by solving Eq. 5 with X = m/N for 
m = m c at which the maximum of K2 is observed. Then, 
the values of K 2 ,K 3 , and K4, calculated from Eq. 5 can 
be verified with those observed at m — m c . 

For all the considered sortings, K 3 = occurs at the 
same position as the minimum of K4, near the maxi- 
mum of the mean cluster size and a point where the 
best power-law fit to the fragment size distribution is 
observed. Therefore, such a point may be used as an 
alternative or complementary indication of the pseudo- 
critical point. The minimum value of K4 is close to — 1 
for all the sorting parameters and system sizes. 

Figs. 10 and 11 show that the crossing points, as it was 



for the sorting by q, are shifted from the "critical" point 
corresponding to p c toward the ordered phase region. For 
Sbound I N they are well localized near the pseudocritical 
point. In this case the pseudocritical point can be ad- 
ditionally characterized by K2 — 0.15 and the position 
Sbound/N ~ 0.84. 

In our simulations events have been generated for uni- 
formly distributed values of the bond probability, p, and 
then grouped in bins of a sorting variable. Using a dif- 
ferent distribution of p changes the spectrum of events in 
a bin. Consequently, quantities such as Ki calculated for 
events within a bin may also change their values. Calcu- 
lations performed for a gaussian distribution of p which 
might simulate experimental conditions have shown that 
this effect is of minor importance. 

IV. CONCLUSIONS 

The largest cluster size distributions have been exam- 
ined within a percolation model on small lattices with 
open boundaries. The dimcnsionless cumulant ratios as 
the normalized variance, K2, the skewness, K3, and the 
kurtosis, K4, of the distribution satisfy with a good ac- 
curacy the finite-size scaling in the critical region. In 
particular, Ki are independent of the system size near 
the critical point (crossing points). This feature has 
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FIG. 12: Parameters of the "critical" point as a function of 
the system size when events are binned according to the mul- 
tiplicity (left column) and Sbound (right column). The points 
show percolation values, the lines represent approximations 
by Eq. (5). 



TABLE I: The coefficients of Eq. (5) for events grouped in 
bins of the multiplicity (upper part) and Sbound (lower part) . 



X 


a 


b 


c 


m/N 


0.352 


28. 


-9.6 


K 2 


0.187 


44. 


6.8 


K 3 


1.03 


19. 


15. 


K 4 


0.91 


26. 


47. 


bound / N 


0.777 


23. 


6. 


K 2 


0.27 


7. 


2.6 


K, 


1.27 


8. 


10. 


K 4 


1.74 


16. 


53. 



been explored in phase transition studies as the cumu- 
lant crossing method, particularly for the kurtosis in the 
form of the Binder cumulant [34| . To our knowledge this 
method has not been applied in analyzing multifragmen- 
tation data. However, it would require a wide range of 
system sizes, which is difficult to realize in nuclear reac- 
tions. Moreover, the presence of crossing points has an 



unambiguous interpretation when events are sorted ac- 
cording to the control parameter. In practice, events are 
grouped with an inevitable dispersion over the control 
parameter. This blurs the scaling behavior; the crossing 
points may appear in a wide range of the sorting param- 
eter away from the "critical" point. These remarks apply 
also for the A = 1 scaling law, since the occurrence of the 
finite-size scaling for the cumulant ratios is a necessary 
condition for this scaling law. The model shows that the 
A-scaling method fails for systems in normal phases. The 
largest cluster size fluctuations in the disordered phase 
cannot be described by a A scaling, whereas the limiting 
A = 1/2 scaling in the ordered phase is violated in small 
systems with open boundaries. 

The percolation transition in such systems can be, 
however, identified by examining some distinct features 
of the finite-size scaling functions of Ki , which are not sig- 
nificantly affected by corrections to scaling in small sys- 
tems and by a dispersion of the control parameter when 
events are sorted according to various measurable quanti- 
ties. The maximum of Ki approximately corresponds to 
the location of the true critical point. The absolute val- 
ues of Ki at this point have been determined for sortings 
by the control parameter, the multiplicity, and Sbound, 
providing complementary characteristics. Coincidentally 
with the maximum of the mean cluster size (pseudocrit- 
ical point) and the power-law behavior of the fragment 
size distribution, one observes K3 = and a minimum 
value of K4 ~ — 1. If the quantity Sbound is used for sort- 
ing events, this point can be additionally characterized 
by Ki ~ 0.15, and its location is related to the system 
size as Sbound — 0.84N. The analysis does not require the 
knowledge nor variation of the system size, which is not 
well controlled in nuclear multifragmentation. It allows 
to estimate the system size at the critical and pseudo- 
critical points. 

It will be interesting to confront these predictions with 
multifragmentation data, in particular with the Aladin 
data, in which the sorting parameter Zbound and the 
charge of the largest fragment are well determined in a 
wide range of the excitation energy. It would be also 
instructive to perform similar analysis with other mod- 
els which are known to contain or not contain a phase 
transition or critical behavior. Using appropriate system 
sizes, boundary conditions and event sortings in model 
simulations is an important requirement. 

This work was supported by the Polish Scientific Re- 
search Committee, Grant No. 2P03B11023. 
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